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Summary 

The SOUSSA surface-panel computer program 
for calculating aerodynamic forces for potential flow 
was applied to a thin fighter-type wing model with 
and without a tip-mounted missile. Steady and un- 
steady lift and pitching-moment results are presented 
for the Mach number M range 0.6 to 0.9 and are 
compared with results from experiments and from 
PANAIR and related computer programs. In gen- 
eral, the SOUSSA program, the experiments, and the 
three PAN AIR (and related) programs give lift and 
pitching-moment results which agree at least fairly 
well, except for the unsteady experimental moment 
at M = 0.8 for the clean wing and the unsteady 
moment at M — 0.8 calculated with a PANAIR- 
predecessor program for the wing with tip-mounted 
missile. 

Introduction 

There has been a continuing quest for more (and 
more detailed) accuracy in determining both steady 
and unsteady aerodynamic forces on wings and 
bodies in the various Mach number ranges. Integral- 
equation methods have been developed which ac- 
count for the actual body shape, thickness, and po- 
sition (orientation) without resorting to in-plane or 
on-axis approximations. A leading approach is to 
evaluate the required integrals by subdividing the 
surface into an array of connected panels, prescribing 
a functional form for the variation of the dependent 
variable (e.g., velocity potential) over each panel, and 
calculating the effect of each panel on all panels. A 
number of such panel methods exist which apply to 
steady flow over a motionless body. 

In contrast, for unsteady motions there are few 
calculation methods available, despite the fact that 
the unsteady capability is of growing importance 
for aeroelastic behavior, control motions, gusts, and 
other transient considerations. The present report 
concerns applications and validation of a calculation 
method which does account for unsteady motions and 
deformations. For inviscid isentropic flow, Morino 
(1974, 1980) employed a generalized Green’s function 
method which transforms the full velocity-potential 
differential equation into an integral-differential de- 
lay equation whose solution for linearized subsonic 
flow is implemented in the prototype computer pro- 
gram SOUSSA Pl.l (Morino 1980 and Smolka et al. 
1980). This program employs a velocity potential 
that is spatially constant over each panel, and it is 
thus termed a “zero-order” method. 

An earlier validation and critique of SOUSSA 
Pl.l was provided by Yates et al. (1982) for five sim- 
ple wing planforms and for one wing-with-fuselage 


configuration. The purpose of the present report is 
to extend the validation to a more geometrically com- 
plex and more realistic configuration, that of the in- 
strumented F-5 fighter wing model both without and 
with a tip-mounted launcher and missile, for which 
the experimental results appear in the reports by 
Tijdeman et al. (1978, 1979). Effects of panel ar- 
rangement and density are studied. 

In order to provide comparisons with other panel 
methods, calculated results from two PAN AIR pro- 
grams and from a predecessor of PANAIR are also 
presented. A PAN AIR program in use at the Langley 
Research Center (see Derbyshire and Sidwell 1982) 
has the capability to calculate only steady forces and 
was used in calculations for the clean wing; these re- 
sults are denoted by PANAIR(S). Another PANAIR 
program, a pilot code (see Dusto and Epton 1980), 
was used to calculate both steady and unsteady com- 
parative results for the clean wing; these results are 
denoted by PANAIR(U). From a higher-order-panel 
method program (predecessor of PANAIR; see Dusto 
1980), steady and unsteady results are included for 
comparison of both the clean wing and the wing 
plus tip-mounted launcher and missile body. These 
three programs solve the linear integral equation of 
unsteady or steady small-disturbance subsonic po- 
tential flow. The wing-body surface is divided into 
panel, and the surface boundary conditions are im- 
posed on the actual steady mean panel locations. 
On each panel the velocity potential is approxi- 
mated by linearly distributed sources and quadrat- 
ically distributed doublets; hence, the term “higher- 
order method” is used. The strengths of sources and 
doublets are solved by satisfying the boundary con- 
ditions. Calculations were made in the Mach number 
range of 0.6 to 0.9. Results are presented in figures 
and listed in tables. 

Symbols 

Subscript Z appears in Tijdeman et al. (1978, 
1979) to identify normal force, which differs only 
insignificantly from lift force in this investigation. 


C L a ’ C *Z a 

lift-curve slope, per radian 

Cm 

pitching-moment coefficient about root 
midchord axis 

/-Y* 

tip-body contribution to pitching- 
moment coefficient about axis 
at strain-gage “balance center” 
(i.e., 75.0- percent root chord) 

Cm o, > C * Ma 

pitching-moment-curve slope about 
root midchord axis, per radian 


Cp pressure coefficient 

Cz lift coefficient 

C* z tip-body contribution to lift coefficient 

c local chord, in. 

Im( ) imaginary part of ( ) 

i imaginary number 

M free-stream Mach number 

Re( ) real part of ( ) 

x local streamwise coordinate, in. 

a angle of attack, deg 

rj spanwise coordinate normalized to 

semispan 


Analysis 

The SOUSSA Program 

For small perturbations and subsonic Mach num- 
bers, the integral equation (2-52) of Morino (1980) 
expresses the relationship between the perturbation 
velocity potential $ and downwash at the surface 
of the body as follows (in the symbols of that report): 
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The velocity potential at body surface point P* at 
present time T* is given on the left-hand side of the 
equation. On the right-hand side are all the effects 
from perturbations that originated in past time T 
from the body and wake surfaces, E^ and Eyy re- 
spectively. The quantity R is the Prandtl-Glauert 
transformed distance traveled by a perturbation in 
reaching P* . The first integral represents a distri- 
bution of sources, the strength of which is the local 
value of downwash on the body surface. The second 
integral represents the effect from a distribution of 
doublets which arise from the interrelation of velocity 


potentials on the body surface, and the third integral 
represents the effect from “wakelets,” the doublets 
distributed on the wake. 

The body surface E# is discretized by quadri- 
lateral panels with hyperboloidal surfaces. The ve- 
locity potential on each panel is assumed to be spa- 
tially constant and of unknown complex amplitude. 
The values of these complex amplitudes are obtained 
from the solution of the set of simultaneous equa- 
tions, one for each panel. The program originally 
provided for a finite-panel wake of arbitrary shape, 
but it was much more economical and satisfactory 
in the present usage to analyze the wake as flat and 
extending downstream to infinity in strips from the 
trailing edge, subject to the Kutta condition, and 
with properties sinusoidally varying downstream ac- 
cording to the past history of harmonic motion. An 
addendum to the SOUSSA manual by Cunningham, 
Desmarais, and Yates (1982) describes the added 
data which control the wake treatment used herein. 

The most general configuration to which any ver- 
sion of SOUSSA has been applied is the Space Shut- 
tle. The full-span vehicle was divided into 764 pan- 
els. Symmetric and antisymmetric subsonic flutter 
analyses were made with 60 vibration modes, includ- 
ing control-surface modes. Also the external tank 
and solid-rocket boosters were modeled and the cal- 
culations were obtained for the launch configura- 
tion. This unpublished work was described by Yates 
(1986). 

A variety of more recent experiences, as well 
as a critique of the program, is given by Yates 
et al. (1982). Five different wing configurations 
and one wing-body configuration were analyzed, and 
steady and unsteady pressure distributions, general- 
ized aerodynamic forces, and flutter-boundary calcu- 
lations were presented. The resulting comparisons 
with experimental results and other analyses were 
mostly very good. 


PANAIR and Predecessor Programs 

The fundamental integral equation for the linear 
steady perturbation velocity potential of PANAIR 
was obtained by applying Green’s third identity to 
convert Laplace’s equation, and it appears as equa- 
tion (3.2.7) of Magnus and Epton (1980), namely 
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where <f>(P ) is the perturbation velocity potential 
at point P on the surface S of the body plus the 
wake. At arbitrary integration points on 5, R is 
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the “compressible” distance to P, h is the unit nor- 
mal, <7 and (i are the unknown source and doublet 
strengths, respectively, and y the gradient opera- 
tor. On the other hand, for the PANAIR-predecessor 
higher-order-panel method used by Dusto (1980), the 
integral equation derives from the application of the 
Prandtl-Glauert transformation to Helmholtz’s equa- 
tion for the unsteady potential function. From pages 
3 and 4 of Dusto and Epton (1980) the integral equa- 
tion is (sic) 
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where Q is the surface integration point, Q = 

< jj/(a/3 2 ) where u is circular frequency, a is speed of 
sound, and /? 2 = (1 — M 2 ), cr* and fi * are the complex 
amplitudes of source and doublet strengths, respec- 
tively, and S', P, and n, have the same meanings as in 
the preceding equation. The close parallelism of the 
two integral equations above is apparent. The differ- 
ence in signs comes from the chosen positive senses 
of the unit normal n. There also exists a PANAIR 
pilot code that includes the unsteady oscillation ef- 
fects, and with this code a few results were calculated 
for comparison. 

The body surface (and optionally the wake) is 
divided into panels, and over each panel there is per- 
mitted a linear variation of source strength and a 
quadratic variation of doublet strength. For each 
point P the surface integrations are subdivided into 
the contributions from each panel per unit value of 
the strengths. The result is a set of simultaneous 
equations whose left-hand side consists of a matrix 
of velocity-potential influence coefficients times a col- 
umn vector of unknown source and doublet strengths. 
The right-hand side of the equation set implements 
the boundary and the Kutta conditions. 

Configuration Definition 

For the present analysis the configuration studied 
is the F-5 fighter wing model for which experimental 
results are reported in the reports by Tijdeman et al. 
(1978, 1979). The half-span wing (see fig. 1) has a 
leading-edge sweep of 31° 55 ; , a taper ratio of 0.308, a 
full-span aspect ratio of 2.98, and a 4.8-percent-thick 
airfoil. The airfoil section has a slight nose-down 
camber over the forward 40 percent of the chord but 
is symmetric over the aft 60 percent. The clean wing 
was tested with a rounded tip fairing. For the tests 
with tip-mounted body the tip fairing was replaced 
with a long, slender missile launcher, and the missile 
body was mounted against the launcher with no gap. 

The clean wing, including the tip fairing, is 
sketched in figure 1. Figure 2 diagrams the wing with 


missile launcher and missile body. Also shown in fig- 
ure 2 is the panelling of the wing and the tip body as 
analyzed with the higher order panelling method and 
reported in Dusto (1980). This panelling is based on 
a nearly cosine distribution plus certain additional 
subdivisions and is termed the “13 X 13 cosine mod- 
ified” panelling. (See tables 1 and 2 of Dusto 1980.) 
Figure 3 shows perspective views of the wing with tip 
body panelled for SOUSSA analysis. All 634 panels, 
including the far-side ones, are shown. The wing it- 
self is divided into panels with a 10 x 10 unform distri- 
bution. Figure 3(a) is an overall view, and figure 3(b) 
is a close-up of the tip body. The body panelling at 
the tip juncture had to be the same streamwise, and 
the rest of the tip body was divided into panels with 
about the same fineness, which was thought to be 
adequate. In general, the 434 tip-body panels were 
approximately square. 


Results and Discussion 

The SOUSSA results were calculated mainly for 
comparison with the available experimental data at 
subsonic Mach numbers of 0.6 and 0.8 and extending 
to M — 0.9 for the wing with the tip body, for which 
M = 0.8 was not tested. Other available panel- 
method results from PAN AIR and PAN AIR- related 
programs are included for comparison. Dusto’s 
(1980) results for the clean wing at M — 0.6 and 
0.8 and for the wing with tip body at M = 0.8 
are included. Also, as stated in the Introduction , a 
few other calculations, designated PANAIR(S) and 
PANAIR(U), are presented for the clean wing. The 
results are organized as indicated in the table on 
page 4. All figures except 7(b) include experimental 
results. All experimental and calculated data pre- 
sented in the figures are also tabulated. Table I for 
the clean wing and table II for the wing with tip body 
list lift and moment and their slopes with respect to 
angle of attack. 

Results are given for the clean wing with tip 
fairing in figures 4 and 5 for M — 0.6 and 0.8, and for 
the wing plus tip-mounted launcher and missile body 
(no fins) in figures 6 and 7 for M = 0.6, 0.8, and 0.9. 
The steady lift and moment coefficients are given for 
a = -0.5°, 0°, 0.5°. The unsteady results are the 
complex amplitudes of lift- and moment-curve slopes 
(per radian) for small-amplitude (±0.11° for clean 
wing, ±0.52° with tip body) pitching oscillations 
about the root midchord axis. For the wing with 
tip body, the lift and moment results are given for 
the entire configuration and separately for the tip 
body itself. The experimental results were obtained 
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on the wing by chordwise integration of pressures for 
the section values and then by spanwise integration 


of the section values for total results. The tip-body 
contributions are from strain-gage data. 


Figure 

Method 

Panelling 

Variable 

(«) 

Motion 

Clean wing 

4(a) 

SOUSSA 

10 x 10, uniform 

) 

Steady 

4(a) 

SOUSSA 

12 x 12, uniform 

p 

Steady 

4(a), 4(b) 

SOUSSA 

13 x 13, cosine, modified 

b 

Steady 

4(a), 4(b) 

Dusto 

13 x 13, cosine, modified 

> M 

Steady 

4(a) 

PANAIR(S) 

13 x 13, cosine, modified 

Vp 

Steady 

4(a) 

PANAIR(S) 

10 x 10, uniform 

) p 

Steady 

4(b) 

PANAIR(U) 

13 x 13, cosine, modified 

> P M 

Steady 

4(b) 

SOUSSA 

13 x 13, uniform 

J a , 1V1 

Steady 

5(a) 

SOUSSA 

10 X 10, uniform 

I P 

Unsteady 

5(a) 

SOUSSA 

13 X 13, cosine, modified 

b 

Unsteady 

5(a), 5(b) 

Dusto 

13 x 13, cosine, modified 

^ M 

Unsteady 

5(a) 

PANAIR(S) 

13 X 13, cosine, modified 

J 

Steady 

5(b) 

SOUSSA 

13 x 13, uniform 

1 \4 

Unsteady 

5(b) 

PANAIR(U) 

13 x 13, cosine, modified 

[ 1V1 

Unsteady 

Wing with tip body 

6(a), 6(b) 

SOUSSA 

10 x 10, uniform (wing) 


Steady 

7(a)-7(c) 

SOUSSA 

10 x 10, uniform (wing) 


Unsteady 

7(b) 

Dusto 

13 x 13, cosine, modified 


Unsteady 


a P Panelling; M — Method. 


Clean Wing 

Steady . Figures 4(a) for M = 0.6 and 4(b) 
for M = 0.8 show steady lift coefficient Cz and 
pitching-moment coefficient Cm (about root mid- 
chord axis) for a — -0.5°, 0°, and 0.5°. Here and for 
other steady results the sloping lines connect selected 
points for the three angles of attack. The experimen- 
tal values (square symbols) are from test runs 136, 
137, and 138 (M = 0.6) and runs 145, 146, and 147 
(M = 0.8) of Tij deman et al. (1978, 1979). It is 
not known why the experimental values vary nonlin- 
early at these small angles of attack. The nonlin- 
earity may be caused by the combination of forward 
camber and the relatively sharp leading edge, which 
results in sizeable down loads over the foremost part 
of the wing, even at a = 0.5°. (See fig. 4(c) discussion 
below.) For M = 0.6 (fig. 4(a)), there are SOUSSA 
results for three different wing panellings: 10 x 10 
uniform (i.e., 10 chordwise and 10 spanwise), 12 x 12 
uniform, and 13 x 13 cosine modified. The latter 
panelling is that employed in Dusto (1980) (as illus- 
trated in his fig. 18 for clean wing with tip fairing) 


and is used for comparison. These three SOUSSA 
results are closely bunched together, thereby indi- 
cating a substantial convergence with respect to the 
panelling. The SOUSSA lift and moment coefficients 
are shown to be a little higher than the experimental 
values (except for a = —0.5°), with the slopes of the 
lift curves being about the same and the slope of the 
SOUSSA moment curves being a little higher than 
the experimental moment curve. Dusto’s results are 
given by the open circles, and both the lifts and the 
moments are slightly above the SOUSSA results. It 
is not known why the DUSTO lift at a = 0° is not 
negative. The steady PAN AIR program was applied 
for the same 13 x 13 panelling used by Dusto. The 
resulting values (open diamond symbols, designated 
PANAIR(S)) fall among the SOUSSA results. How- 
ever, the PANAIR(S) code with the 10 x 10 uniform 
panelling results in both lift and moment coefficients 
which are sharply different from the other results. 
This difference is unexpected and is thought to be 
caused by the relatively large leading-edge panels and 
the associated calculation of lifting pressure behind 
the leading-edge peak. 
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Figure 4(b) compares calculations from three dif- 
ferent codes with experimental values for M — 0.8. 
The Dusto (1980) results with the 13 x 13 cosine 
modified panelling and the SOUSSA results with the 
13 x 13 uniform panelling have the same relative posi- 
tions as for M — 0.6; that is, the analytical results are 
moderately above the experimental values, except at 
a = 0.5°, and the Dusto results are near and above 
the SOUSSA results. The unsteady PANAIR(U) pi- 
lot code results are given by the solid circles. They 
are very close to the SOUSSA lift results but are dis- 
tinctly below the experimental moment results. All 
the calculations overpredicted Cz at a — 0° and its 
slope somewhat, as is typical for potential flow 
methods. 

Figure 4(c) gives a comparison of SOUSSA and 
experimental pressure distributions at three span sta- 
tions for M = 0.8 and a = 0.5°. The SOUSSA 
results from a 13 x 13 uniform panelling are plot- 
ted as straight-line segments connecting the panel- 
center values. At the inboard and midspan sta- 
tions, SOUSSA pressures are mostly below the ex- 
perimental values over the forward part of the chord, 
although lifting pressure differences are well pre- 
dicted. Outboard, however, forward lifting pressure 
arc slightly underpredicted. As mentioned previ- 
ously, both calculation and experiment indicate size- 
able down loads on the foremost part of the wing all 
across the span, even for a — 0.5°. 

Unsteady . Figures 5(a) for M — 0.6 and 5(b) for 
M — 0.8 show the slopes (per radian) of the complex 
amplitude of the lift coefficient and of the pitching- 
moment coefficient for small amplitude (±0.11°) 
pitching oscillations about the zero mean angle of at- 
tack. In part II of their report, Tijdeman et al. (1979) 
list the oscillation amplitude as 0.11°. For both Mach 
numbers, in addition to the quasisteady values there 
are experimental values for two pitching frequencies, 
20 and 40 Hz, a Dusto (1980) result for 40 Hz, and 
SOUSSA results for 0, 20, and 40 Hz. Here and for 
other unsteady results the lines from the origin to se- 
lected points serve to display the phase angle of that 
complex result. The reduced frequency for 40 Hz is 
0.399 for M — 0.6 and is 0.307 for M = 0.8, based on 
root half-chord. Each calculation used the matching 
experimental reduced frequencies. The quasisteady 
results are inferred from the slope of the static re- 
sults. For M — 0.8, the PANAIR(U) results also are 
presented. 

For M — 0.6 (fig. 5(a)), the calculated and ex- 
perimental results display fair to good agreement. 
In comparison to experimental results, calculated lift 
and moment magnitudes are about 25 percent larger, 
moment phase angles agree well, SOUSSA lift phase 


angles are a few percent less (about 5°) for 20 Hz 
but are 20 to 30 percent less at 40 Hz, and the 
Dusto (1980) lift phase angle is about 45 percent 
(about 12°) less at 40 Hz. The agreement between 
the SOUSSA results and the experimental results im- 
proves significantly as the panelling is made finer, 
thereby indicating that the convergence with respect 
to panelling is not as well established for the unsteady 
case as for the steady case. 

Figure 5(b) gives the calculated and experimen- 
tal results for M = 0.8. Compared with the exper- 
imental lift-slope magnitudes, SOUSSA magnitudes 
are about 30 percent larger, PANAIR(U) magnitudes 
are about 25 percent larger, and Dusto (1980) magni- 
tudes are 35 to 45 percent larger. All the calculated 
lift-slope phase angles are smaller than the experi- 
mental angles — the SOUSSA by 40 percent (about 
9° at 40 Hz), the PANAIR(U) by 55 percent (about 
12°), and the Dusto (1980) by 65 percent (about 15° 
at 40 Hz). For the moment slope the experimen- 
tal quasisteady moment slope is positive (i.e., in the 
pitch-divergent sense) as for all other cases studied, 
but for 20 and 40 Hz the moment slope has shifted 
over into the third quadrant, indicating a real part 
that opposes pitch displacement. This shift from the 
0-Hz to the 20-Hz moment slope represents a rear- 
ward shift of center of lift by 8 percent of root chord. 
The negative imaginary part continues to damp the 
pitching velocity. 

In contrast to the experimental results, all the 
calculated moment slopes are in the fourth quad- 
rant. For all frequencies the SOUSSA moment mag- 
nitudes are twice the experimental values, and the 
PANAIR(U) magnitudes are even higher. In com- 
parison with SOUSSA, the Dusto (1980) result for 
40 Hz is very close. It is notable and somewhat sur- 
prising that the PANAIR(U) code, with its higher or- 
der panels, has not given a plainly better result than 
the SOUSSA code with its “zero order” (i.e., con- 
stant potential) panels. Here again, as in figures 4(a) 
and 4(b), all the calculations substantially overpre- 
dict the magnitudes of both C^ a and Cj^ a at both 
Mach numbers. 

Figure 5(c) shows distributions of lifting-pressure 
slope results from SOUSSA and experiment at three 
span stations at M — 0.8 for pitching motion about 
a = 0°. The real and imaginary parts of the 
complex amplitude (per radian) are plotted. The 
SOUSSA results are plotted as straight-line segments 
connecting the 13 panel-center values. The chordwise 
distribution of experimental lifting-pressure slopes 
show a strong reversal at the most forward data point 
for r? = 0.18, a less strong reversal for r] = 0.51, 
and no reversal for t] = 0.87. The SOUSSA results 
show no reversals with the relatively coarse panelling 
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used near the leading edge. Aft of the most forward 
data point the SOUSSA results agree well for the 
imaginary part and somewhat less well for the real 
part, with some underprediction of the real part for 
x/c = 0.20 to 0.40, and some over prediction aft for 
x/c = 0.70 to 0.90. 

Wing With Tip Body 

The remaining results are for the wing plus tip- 
mounted launcher and missile body. For such a ge- 
ometrically complex configuration the SOUSSA pro- 
gram requires a large amount of tightly organized and 
coordinated geometrical input data, which is gener- 
ated by a preprocessor program. The preprocessor 
program used was not efficiently organized, and as 
a result only a limited number of panels could be 
generated within the available machine core memory 
(131072 words in a CDC 6600 computer). Because 
of this preprocessor limit, the wing plus tip body 
was divided into 634 panels, with the wing upper 
and lower surfaces each containing 10 panels chord- 
wise by 10 spanwise. Relatively few additional panels 
could have been accommodated, and hence the fins 
could not be added. 

In each figure the total configuration results ap- 
pear at the left, and the contributions from the tip 
body appear at the right. The reference wing area 
is the clean-wing area for all coefficients. The mo- 
ment reference length is the mean geometric chord, 
0.4183 m. The tip-body pitching moment C* M is ref- 
erenced about an axis at the strain-gage balance cen- 
ter at 75.0 percent of the root chord. 

Steady . The steady results are given in fig- 

ures 6(a) for M — 0.6 and 6(b) for M = 0.9. 
The comparison between the experimental and the 
SOUSSA results is good. In fact, the comparison is 
mostly better than for the clean wing in figure 4. The 
agreement for the tip-body contributions is especially, 
good. For M = 0.6, the SOUSSA lift magnitudes at 
a = ±0.5° and the lift-curve slope are increased by 
about 10 percent over that of the clean wing, while 
the lift on the tip body is just under 3 percent of the 
total. This indicates a significant “end plate” effect 
of the tip body acting on the wing. The upsetting 
moment is decreased by nearly 10 percent because 
the increased lift is added largely behind the pitch 
axis, and it contributes in the restoring sense (i.e., 
negatively). For M = 0.9 the comparisons are very 
similar to those for M = 0.6. 

Unsteady . Figures 7(a) for M ~ 0.6, 7(b) for 
M = 0.8, and 7(c) for M = 0.9 show the unsteady 
results for the total configuration on the left and the 


tip-body contributions on the right. For M = 0.6, 
the magnitudes of the total-configuration lift-curve 
slopes from SOUSSA for all three frequencies are 
about 10 percent greater than for the clean wing 
of figure 5(a). The tip-body contributes less than 
3 percent of this lift, thereby again indicating an 
end-plate effect of the tip body. This increased lift 
near the tip acts to resist the pitch displacement 
and causes the SOUSSA total-configuration moment- 
curve slope to be about 10 percent less than for the 
clean wing. The magnitudes of the SOUSSA total- 
configuration lift curve, total-configuration moment 
curve, and tip-body moment curve are about 25 per- 
cent greater than the 20-Hz experimental values, and 
the SOUSSA tip-body lift curve is slightly less than 
the experimental. 

For M = 0.8 (fig. 7(b)), there are no experimen- 
tal results but there are comparisons of the quasi- 
steady and 40-Hz results from SOUSSA and Dusto 
(1980). The comparisons are fairly good except for 
the pitching moment. The Dusto (1980) pitching- 
moment slope in the third quadrant is from his ta- 
ble 9 for configuration C4. The reason for the dif- 
ferent quadrant is not known, although Dusto states 
that the Kutta condition was not achieved in the un- 
steady case, thereby implying that an improvement 
was needed. For the other three comparisons, the 
SOUSSA phase is consistently greater by a consider- 
able percentage margin, although all the phase angles 
are rather small. 

For M — 0.9 (fig. 7(c)), there are comparisons of 
the quasisteady and the 20-Hz SOUSSA and exper- 
imental results. At 20 Hz the reduced frequency is 
0.138 (reference length is the root half-chord). The 
pitch amplitude of the experiment is 0.53°. The 
SOUSSA lift-curve and moment-curve magnitudes 
are somewhat greater than experimental magnitudes 
for the total configuration but are somewhat less 
for the tip-body contribution. For all results except 
the total-configuration moment curve, the SOUSSA 
phase angle is about double that of the experiment, 
although again all phase angles are small. For the 
total-configuration moment curve, the phase angle 
agreement is good. Notably, the SOUSSA pitching 
moment is in the fourth quandrant of the complex 
plane for all three Mach numbers. 

To summarize briefly, in general the three dif- 
ferent analyses and the experiment give lift and 
pitching-moment results which agree among them- 
selves at least fairly well, except for the unsteady 
experimental moment on the clean wing at M = 0.8 
and the unsteady Dusto-calculated moment on the 
wing plus tip body at M — 0.8. SOUSSA, with its 
zero-order panel approximation (i.e., uniform veloc- 
ity potential on each panel), has performed rather 
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well in comparison with experimental results and 
with PANAIR, a higher order panel method. 

Concluding Remarks 

Application of a generalized Green’s function 
method to the velocity-potential differential equa- 
tion transforms it into an integral-differential delay 
equation which can be implemented for computation 
by surface panelling. The prototype surface-panel 
computer program SOUSSA Pl.l, which is based on 
this concept, calculates subsonic steady and unsteady 
pressure distributions and integrated forces on bodies 
of arbitrary shape and deformation. 

In an earlier validation and critique of the 
SOUSSA panel method, calculations were made for 
six simple configurations. The purpose of the present 
report was to extend the validation to a geometrically 
more complex and more realistic configuration, that 
of an instrumented F-5 wing model with and without 
a tip-mounted missile. 

Calculated lift and moment results were com- 
pared with experimental results. Steady and un- 
steady results were obtained at Mach numbers M 
of 0.6, 0.8, and 0.9 at steady pitch angles of 0° and 
TO. 5° and for small pitching oscillations at 20 and 
40 Hz. For the clean wing, where variations of the 
SOUSSA panelling were used, very good convergence 
was indicated at M — 0.6 for steady flow, but for un- 
steady pitching oscillations the modest differences at 
40 Hz indicated less assured convergence with respect 
to panelling. For the wing with tip body, panelling 
variations were not performed. 

The agreement of SOUSSA calculations with ex- 
perimental results ranged from very good to fairly 
good with the exception of the unsteady pitching- 
moment slope on the clean wing at M — 0.8, for 
which the experiment may conceivably be in error. 
The SOUSSA results also agreed reasonably well 
with most of the results from three different PANAIR 
and PANAIR-related programs. A steady-flow-only 
PANAIR program produced distinctly poorer agree- 
ment for lift and pitching moment for uniform chord- 
wise panelling (i.e., with a “large” leading-edge 
panel) than for a modified cosine panelling distribu- 
tion (with a much smaller leading-edge panel). 

In general the SOUSSA program, the experi- 
ments, and the three PANAIR-related programs gave 
lift and pitching-moment results which agreed at 
least fairly well, except for the unsteady clean- wing 
experimental moment at M — 0.8 and a moment cal- 
culated with a PANAIR predecessor at M — 0.8 on 
the wing with tip body. SOUSSA, with its zero-order 
panel approximation (i.e., uniform velocity poten- 
tial on each panel), has performed rather well com- 


pared with experiment and with a higher-order panel 
method, PANAIR. 


NASA Langley Research Center 
Hampton, VA 23665-5225 
June 12, 1987 
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Table I. Results for Clean Wing 
(a) M = 0.6, steady (fig. 4(a)) 


Source 

a, deg 

Cz 

Cu 

Cm 

Cm* 

Experiment, run 136 

-0.5 

-0.0342 


-0.0116 


Experiment, run 137 

0 

-.0093 

3.06 

-.0096 

0.241 

Experiment, run 138 

.5 

.0192 


-.0074 


SOUSSA, 10 x 10, 

-0.5 

-0.0367 


-0.0116 


uniform 

0 

-.0050 

3.632 

-.0088 

0.323 


.5 

.0267 


-.0060 


SOUSSA, 12 x 12, 

-0.5 

-0.0376 


-0.0113 


uniform 

0 

-.0059 

3.627 

-.0085 

0.320 


.5 

.0257 


-.0057 


SOUSSA, 13 x 13, 

-0.5 

-0.0348 


-0.0112 


cosine, modified 

0 

-.0031 

3.621 

-.0084 

0.324 


.5 

.0285 


-.0055 


Dusto (1980, table 8) 

-0.5 

-0.0318 


-0.00935 



0 

.0010 

3.762 

-.00721 

0.245 


.5 

.0338 


-.00508 


PANAIR(S), 

-0.5 

-0.0119 


-0.0164 


10 x 10, uniform 

0 

.0198 

3.64 

-.0145 

0.218 


.5 

.0516 


-.0126 


PANAIR(S), 13 x 13, 

—0.5 

-0.0343 


-0.0106 


cosine, modified 

0 

-.0024 

3.66 

-.0087 

0.227 


.5 

.0296 


-.0067 



(b) M = 0.8, steady (fig. 4(b)) 


Source 

a, deg 

Cz 

Cl* 

Cm 

Cm* 

Experiment, run 145 

-0.5 

-0.0318 


-0.0125 


Experiment, run 146 

0 

-.0082 

3.09 

-.0117 

0.195 

Experiment, run 147 

.5 

.0222 


-.0091 


SOUSSA, 13 x 13, 

-0.5 

-0.0390 


-0.0136 


uniform 

0 

-.0040 

4.013 

-.0105 

0.360 


.5 

.0311 


-.0073 


Dusto (1980) 

-0.5 

-0.0346 


-0.0109 



0 

.0014 

4.13 

-.0087 

0.260 


.5 

.0375 


-.0064 


PANAIR(U), 13 x 13, 

-0.5 

-0.0403 


-0.0175 


cosine, modified 

0 

-.0063 

3.88 

-.0142 

0.390 

j 

.5 

.0275 


-.0107 
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Table I. Concluded 
(c) M = 0.6, unsteady (fig. 5(a)) 


Source 

Frequency, Hz 

C L a 

Cm* 

Experiment 

0 

3.060 + 0.0? 

0.241+0.0? 

Experiment, run 382 

20 

3.097 + 0.702? 

.280-0.195? 

Experiment, run 383 

40 

3.075 + 1.621? 

.294 - 0.372? 

SOUSSA, 10 x 10, 

0 

3.632 + 0.0* 

0.323 + 0.0? 

uniform 

40 

3.514 + 1.275* 

.352 - 0.449? 

SOUSSA, 13 x 13, 

0 

3.621 + 0.0? 

0.324 + o.o?; 

cosine, modified 

20 

3.568 + 0.680? 

.327 - 0.253? 


40 

3.540 + 1.457* 

.353 - 0.502? 

Dusto (1980), 13 x 13, 

0 

3.76 + 0.0? 

0.245 + 0.0? 

cosine, modified 

40 

3.78 + 1.084? 

.354 + 0.464? 

PANAIR(S),13 x 13, 
cosine, modified 

0 

3.661 + 0.0? 

0.227 + 0.6? 


(d) M = 0.8, unsteady (fig. 5(b)) 


Source 

Frequency, Hz 

c L a 

C M a 

Experiment 

0 

3.090 + 0.0? 

0.195 + 0.0? 

Experiment, run 367 

20 

2.938 + 0.609? 

-.039-0.141? 

Experiment, run 368 

40 

3.132 + 1.281? 

-.027 - 0.323? 

SOUSSA, 13 x 13, 

0 

4.013 + 0.0? 

0.360 + 0.0? 

uniform 

20 

3.994 + 0.428? 

.355 -0.281? 


40 

4.044 + 0.922? 

.347 - 0.569? 

Dusto (1980), 13 x 13, 

0 

4.130 + 0.0? 

0.260 + 0.0? 

cosine, modified 

40 

4.483 + 0.531? 

.335 - 0.577? 

PANAIR(U), 13 x 13, 

0 

3.88 + 0.0? 

0.390 + 0.0? 

cosine, modified 

20 

3.85 + 0.34? 

.447 - 0.275? 


40 

3.85 + 0.726? 

.425 - 0.527? 
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Table II. Wing With Tip Body 
(a) M — 0.6, steady (fig. 6(a)) 


Source 

a, deg 

Cz 

C La 

c z 

C L 

L/ql 

Cm 

CMoc 

Cm 

Cm a 

Experiment, run 255 

-0.5 

-0.0471 


-0.00113 


-0.0120 


-0.00075 


Experiment, run 256 

0 

-.0086 

3.97 

-.00002 

0.128 

-.0101 

0.241 

-.00018 

0.0493 

Experiment, run 257 

.5 

.0221 


.00110 


-.0078 


.00039 


SOUSSA, 10 x 10, 

-0.5 

-0.0403 


-0.00119 


-0.0113 


-0.00075 


uniform 

0 

-.0056 

3.98 

-.00020 

0.128 

-.0088 

0.298 

-.00018 

0.0493 


.5 

.0292 


.00080 


-.0062 


.00039 



(b) M = 0.9, steady (fig. 6(b)) 


Source 

a, deg 

Cz 

C L a 


ES 




C*M a 

Experiment, run 249 

-0.5 

-0.0528 


-0.00149 


-0.0171 


-0.00071 


Experiment, run 251 

0 

-.0091 

4.73 



-.0138 



0.058 

Experiment, run 252 

.5 

.0297 




-.0121 


.00030 


SOUSSA, 10 x 10, 

-0.5 

-0.0469 




-0.0154 


-0.00092 


uniform 

0 


4.89 ; 





-.00023 

0.0779 


.5 







.00044 



(c) M = 0.6, unsteady (fig. 7(a)) 


Source 

Frequency, 

TIz 

c u 

C Za 

°M a 

c* 

M a 

Experiment 

0 

3.97 + 0.0* 


0.241 + 0.0 i 

0.0493 + 0.0* 

Experiment, run 258 

20 


.131 + 0.0245* 

.212 - 0.198* 

.0474 + 0.0019* 

SOUSSA, 10 x 10, 

0 

3.977 + 0.0 i 

0.114 + 0.0*' 

0.298 + 0.0*' 

0.065 + 0.0* 

uniform 

20 

3.893 + 0.592* 

.112 + 0.0281 

.301 - 0.245* 

.064 - 0.00135* 


40 

3.824 + 1.321* 

.109 + 0.059z 

.328 - 0.484* 

.063 0.00039* 


(d) M = 0.8, unsteady (fig. 7(b)) 


Source 


C L a 



C Ma 

SOUSSA, 10 x 10, 

B91 

4.483 + 0.0* 


0.323 + 0.0* 


uniform 

- 1 

4.439 + 0.924* 

,129 + 0.042* 

.305 - 0.595* 


Dusto (1980), 13 X 13, 

0 

4.595 + 0.0* 


0.339 + 0.0* 

0.0447 + 0.0* 

cosine, modified 


4.421 + 0.656* 

.165 + 0.0270*' 

-.247 - 0.640* 

.0549 - 0.0083*' 


(e) M = 0.9, unsteady (fig. 7(c)) 


Source 

Frequency, 

Hz 

C L a 


C M a 

C* 

Jvfcjt 

Experiment 

0 

4.73 + 0.0* 

0.164 + 0.0*' 

0.286 + 0.0*' 

0.0579 + 0.0*' 

Experiment, run 253 

20 

4.13 + 0.138* 

.170 + 0.0060*' 

.231 - 0.381* 

.0550 + 0.0075* 

SOUSSA, 10 x 10, 

0 

4.89 + 0.0* 

0.144 + 0.0* 

0.344 + 0.0* 

0.0456 + 0.0* 

uniform 

20 

4.91 + 0.397* 

.145 + 0.0184* 

.306 - 0.449* 

.0442 - 0.0130* 
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(a) Wing with tip body. 

Figure 3. SOUSSA 10 x 10 uniform wing panelling (total of 634 panels). 
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(b) Close-up of tip body. 
Figure 3. Concluded. 
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z 



EXP., RUNS 136.137,138 
PANAIR(S). 10X10 UNIFORM 
PANAIR(S) \ 

Dusto f 13X13 COSINE. MOD. 
SOUSSA ) 

SOUSSA, 12X12 UNIFORM 
SOUSSA, 10X10 UNIFORM 



a, deg 


a, deg 


(a) M = 0.6, lift and moment. 

Figure 4. Steady lift, pitching-moment, and pressure coefficients for clean wing. 
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a, deg 


a, deg 


(b) M = 0.8, lift and moment. 
Figure 4. Continued. 
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□ EXP., RUNS 382,383 
O Dusto. 13X13 COS. MOD. 

SOUSSA, 13X13 COS. MOD. 

& SOUSSA. 10X10 UNIFORM 
Q P AN AIR(S), 13X13 COS. MOD. 



(a) M — 0.6, lift and pitching-moment slopes. 

Figure 5. Unsteady lift, pitching-moment, and pressure slopes for clean wing pitching about a = 0°. 
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© — ^ 0 Hz 



□ EXP., RUNS 367,368 
O Dusto, 1 3X1 3 COS. MOD. 

^ SOUSSA, 1 3X1 3 UNIFORM 
^ PANAIR(U), 1 3X1 3 COSINE, MOD 


] 40 Hz 


□40 



(b) M = 0.8, lift and pitching-moment slopes. 
Figure 5. Continued. 




□ EXP., RUNS 255,256,257 
A SOUSSA, 10X10 UNIFORM 
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(a) M = 0.6. 

Figure 6. Steady lift and pitching-moment coefficients on wing with tip body and contributions from tip body 
itself. 
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□ EXP., RUNS 249,251,252 
A SOUSSA, 10X10 UNIFORM 
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(b) M = 0.9. 
Figure 6. Concluded. 
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Ot 


□ EXP., RUNS 258,259 
A SOUSSA, 10X10 UNIFORM 


lm(C L ) 

a 




(a) M = 0.6. 

Figure 7. Unsteady lift and pitching-moment coefficients on wing with tip body and contributions from tip 
body itself. 
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Re(C M( ) R«( C M a *) 



& SOUSSA, 10X10 UNIFORM 
O Dusto, 13X13 COSINE. MOD 



(b) M = 0.8. 
Figure 7. Continued. 
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PD 


EXP. 

SOUSSA, 10X10 UNIFORM 



(c) M = 0.9. 
Figure 7. Concluded. 
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